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ABSTRACT 

We present optical R-band light curves of five SDSS double QSOs (SDSS J0903+5028, SDSS J1001+5027, SDSS J1206+4332, SDSS 
J1353+1138, SDSS J1335+0118) obtained from monitoring at the Nordic Optical Telescope (NOT) between September 2005 and 
September 2007. We also present analytical and pixelated modeling of the observed systems. For SDSS J1206+4332, we measured 
the time delay to be Ar = 1 16^ days, which, for a Singular Isothermal Ellipsoid model, corresponds to a Hubble constant of 73*? 
km s" 1 Mpc -1 . Simultaneous pixeleted modeling of five other systems for which a time delay has now been previously measured at 
the NOT leads to Hq = 61.5** km s -1 Mpc - '. Finally, by comparing lightcurves of the two images of each system, suitably shifted by 
the predicted or observed time-delays, we found no evidence for microlensing variability over the course of the monitoring period. 
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1. Introduction 

A gravitationally lensed quasar is a quasar (QSO) which is 
lensed by a massive foreground object such as a galaxy or a 
group of galaxies. If the lens is close enough to the line of sight 
then the quasar is strongly lensed and will be multiply imaged. 
Due to the different travel times for each light path, any intrin- 
sic variation of the quasar is observed in the images at different 
times. This time difference, referred to as the time delay, can be 
measured by comparing the light curves of the images. As first 
shown by Refsdal (1964) the Hubble constant can be determined 
from the time delay provided that the mass distribution is known. 
Conversely, for a fixed Hubble constant, the mass distribution of 
the lensing galaxy can be constrained. 

Determining the Hubble constant from the time delays be- 
tween light variations in multiple images of gravitationally 
lensed QSOs is a classical application of lensing in cosmology, 
but despite 25 years of effort, results remain inconclusive. This 
is partly due to incomplete knowledge of the mass distribution 
along the light path to the QSO and partly due to the perturbing 
effects of microlensing on sparsely sampled light curves of the 
underlying long-term variations. 

To address these issues, we have conducted a monitoring 
program at the Nordic Optical Telescope (NOT) with the aim of 
obtaining densely sampled light curves of five lensing systems 
(SDSS J0903+5028, SDSS J1001+5027, SDSS J1206+4332, 
SDSS J1353+1 138, SDSS J1335+01 18). Dense sampling is re- 
quired to quantify the effects of microlensing (Paczyriski 1986; 
Schild 1996; Paraficz et al. 2006) and to maximize chances of 
determining a time delay from a relatively short monitoring cam- 
paign. 

Estimating the Hubble constant using time delay measure- 
ments is strongly dependent on the underlying mass distribution 
and hence the choice of lens model (Oguri 2007). Two differ- 



quasars: individual: SDSS J0903+5028, SDSS J1001+5027, 



ent approaches to modeling lenses are commonly used. The first 
one, the non-parametric method (Saha & Williams 1997), gener- 
ates a large number of models which perfectly fit the data, each 
of them gives a different time delay which can be then averaged. 
For the second method, the analytical method(Keeton 2001), one 
assumes physical properties of the mass distribution of the lens. 
Comparison of the two approaches gives a useful indication on 
the systematic errors in e.g. the Hubble constant determined this 
way. 

Luckily, it seems that simple lens models are very good 
first approximations to the real mass distributions of lenses 
(Koopmans et al. 2006). Therefore we have chosen to use the 
singular isothermal potential to analytically model all the lenses 
we observe. 

In this paper we present the results of a monitoring campaign 
at the NOT of 5 doubly lensed quasars. We have measured the 
time delay of one of the lensing systems, SDSS J 1206+4332, 
by analyzing light curves of the two quasar images obtained 
from six months of monitoring, demonstrating the feasibility 
of short-term monitoring for time-delay measurement. Based on 
our measurement of the time delay and 5 other time delay mea- 
surements previously obtained at the NOT we have estimated the 
Hubble constant. 

Section 2 describes the details of our monitoring campaign 
and section 3 introduces the observed targets. The photometric 
technique based on image deconvolution is described in section 
4 and the light curves of all 5 targets are presented in section 
5. In section 6 the time delay of SDSS J1206+4332 is deter- 
mined. In section 7 we perform a microlensing search in the 5 
quasar light curves and estimate upper limits to the microlensing 
signal. In section 8 we perform analytical and pixeleted model- 
ing. Section 9 is devoted to simultaneous modeling of 5 NOT- 
determined time delay systems with the aim of a joint Hubble 
constant estimate. In section 10 we discuss the results. 
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In the paper we use a fiat ACDM Universe, Q, m = 0.3 and 
O a = 0.7. 

2. Observation 

The observations were gathered from monitoring programs car- 
ried out in the periods September-March 2005/2006, October- 
March 2006/2007 and April-September 2007 at the Nordic 
Optical Telescope (NOT); a 2.5-m telescope located at Roque de 
los Muchachos in La Palma, Spain. The advantage of this tele- 
scope is its fairly flexible scheduling, which made almost nightly 
monitoring possible. Our targets were observed every night un- 
der all three operative modes at the NOT: observer, service and 
technical. However there were severe obstacles preventing fre- 
quent sampling: bad weather, Guaranteed Time programs at the 
NOT and sensitivity of the deconvolution software to imperfect- 
ness of the data. 

The detectors used in the monitoring were chosen in or- 
der to obtain the most frequent sampling and the best image 
quality. Thus, we used ALFOSC (The Andalucia Faint Object 
Spectrograph and Camera; pixel scale 07189) whenever it was 
mounted on the telescope and StanCam (the Stand-by CCD 
Camera, which is permanently mounted at the NOT; pixel scale 
07176), otherwise. The seeing varied from 074 to 370 with 170 
being the most frequent value. 

We observed the objects in the R band, only. The pilot phase 
of the program involved monitoring of a gravitationally lensed 
system SDSS J0903+5028 at z = 3.6 (Johnston et al. 2003) for 
10 minutes every night. This system contains two quasar images 
(see Table 3.1) separated by 278 and aligned on opposite sides 
of the lensing galaxy at the redshift z = 0.388 (Johnston et al. 
2003). 

After the first observing season we decided to observe 
three lensed systems alternately, SDSS J1001+5027 and 
SDSS J1206+4332 together in one night for 5 minutes each, and 
SDSS J0903+5028 every second night for 10 minutes. The third 
observing season was again divided in two groups. We contin- 
ued observing SDSS J1001+5027 and SDSS J1206+4332 for 5 
minutes each every second night, alternately with two new tar- 
gets, SDSS J1335+01 18 (5 min) and SDSS J1353+1 138 (3 min). 
The exposure times were chosen so that the signal-to-noise ra- 
tio of the fainter component of each system would be above 10. 
Finding charts of all the targets are presented in Figures 1 and 2. 

One of the major challenges in monitoring gravitationally 
lensed quasars is the lack of prior knowledge of time delay. Thus, 
time-delay measurements require well-sampled light curves with 
accurate photometry over a period of time substantially longer 
than the time delay. Because the time delay is not known, one 
has to calculate a theoretical time delay before planning the ob- 
servations. 

Assuming that the theoretical prediction of the time delay for 
a given system is a good approximation (which is not always the 
case) there is still a question as to whether the quasar will vary 
during the period of monitoring and what will be the timescale 
and amplitude of the variations. Quasar brightness might vary 
from a day up to years independently of the mass of the black 
hole (Wold et al. 2007). We know that the rapid variability im- 
plies a light source at very small distance from the black hole 
(Webb & Malkan 2000), while variations on long time-scales 
are related to morphological changes of jets on parsec scales or 
to accretion-disk instabilities (de Vries et al. 2006). However, it 
is impossible to predict whether a quasar will vary or not within 
a given period and timescale, which is why the success of each 
quasar monitoring is uncertain. 



Another difficulty with gravitationally lensed quasar ob- 
servations is microlensing by stars in the lensing galaxy. 
Microlensing can change the results by about 0.5 - 1.0 mag- 
nitudes and it is completely unpredictable (Chang & Refsdal 
1979). 

Thus, in order to get a time delay one has to monitor quasars 
with high sampling, so that events like microlensing can be ex- 
tracted from the quasar variability in further analysis. Our mon- 
itoring program was designed to minimize this problem. 

Observation planning, monitoring supervision and image re- 
duction was made by the first author during her stay at the NOT. 
Data reduction was performed using IRAF reduction utilities. 

3. Targets 

The individual targets observed are introduced below. The as- 
trometry and redshifts of the quasar images and the lenses are 
summarized in Table 3.1. 

3.1. SDSS J0903+5028 

SDSS J0903+5028, a doubly lensed quasar system, was discov- 
ered from the Sloan Digital Sky Survey by Johnston et al. (2003). 
Using the ARC 3.5 meter telescope it was found that the sys- 
tem has two quasar images separated by 278 with the lens, a red 
galaxy (z = 0.388), in between. Spectroscopic follow up obser- 
vation at the Keck II telescope proved that the two objects are the 
images of one quasar at z — 3.6. Johnston et al. (2003) concluded 
that other galaxies in the vicinity of the lensing galaxy might be 
gravitationally bound with the lens, adding external shear to the 
lensed system. 

3.2. SDSS J1 001+5027 & SDSS J1 206+4332 

SDSS J1001+5027 and SDSS J1206+4332 were recognized 
as lensing candidates in a strong lensing survey using the 
Sloan Digital Sky Survey. Oguri et al. (2005) photometrically 
and spectroscopically confirmed using the ARC 3.5 meter and 
University of Hawaii 2.2-meter telescopes, that those two tar- 
gets are indeed lensing systems. They reported that SDSS 
J1001+5027 and SDSS J1206+4332 have image separation of 
2786 and 2790 respectively with the source objects being quasars 
atz= 1.838 and z = 1.789. 

Oguri et al. (2005) noted that the two lensing galaxies of 
SDSS J1001+5027 have colors consistent with those of early- 
type galaxies at 0.2 < z < 0.5. SDSS J1206+4332 appears to 
have three lensing galaxies where the main lens, Gl, has an as- 
sociated absorber at redshift z = 0.748, the second galaxy, G2, is 
identifed as a high redshift galaxy z < 0.7 and, the third galaxy, 
G3, is a blue galaxy. Oguri et al. (2005) using the lensmodel 
software (Keeton 2001) modeled the two systems showing that 
they are both strongly affected by the potential of more than one 
galaxy. 

3.3. SDSS J1 335+01 18 

SDSS J1353+1138 is a doubly lensed system, discovered in 
the Sloan Digital Survey by Oguri et al. (2004). The photo- 
metric follow up made at the Subaru 8.2-m and Keck I tele- 
scopes confirmed that the system consists of two gravitationally 
lensed images separated by 1756 with a single lensing galaxy 
in the centre. Spectroscopic observations made at the ESO New 
Technology Telescope (NTT) showed that the A and B compo- 
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Fig. 1. Finding chart of the three systems: Top: SDSS J0903+5028, 
Middle: SDSS J1001+5027, Bottom: SDSS J1206+4332. The refer- 
ence star (SI) and the star used to model the PSF are indicated. 



nents of the system are images of a quasar at redshift z — 1.57. 
Eigenbrod et al. (2006b) reported that the lensing galaxy is a 
low-redshift galaxy with z = 0.44. 



SDSS J0903+5028 


Object 


R.A. (J2000.0) 


Dec. (J2000.0) 




A 


09 03 35.13 


+50 28 20.21 


3.6 


B 


09 03 34.88 


+50 28 18.75 


3.6 


G 


09 03 34.93 


+50 28 19.53 


0.388 


SDSS J1001+5027 


Object 


R.A. (J2000.0) 


Dec. (J2000.0) 


z 


A 


10 01 28.61 


+50 27 56.9 


1.838 


B 


10 01 28.35 


+50 27 58.5 


1.838 




x[arcsec] 


y[arcsec] 


z 


Gl 


1.779 ±0.049 


0.857 ±0.123 


0.2 <z< 0.5 


G2 


1.795 ±0.088 


-0.700 ± 0.053 


0.2 < z < 0.5 


SDSS J1206+4332 


Object 


R.A. (J2000.0) 


Dec. (J2000.0) 


z 


A 


12 06 29.65 


+43 32 17.6 


1.789 


B 


12 06 29.65 


+43 32 20.6 


1.789 




x[arcsec] 


y[arcsec] 




Gl 


-0.664 ±0.137 


1.748 ±0.028 


0.748 


G2 


1.320 ±0.147 


5.999 ±0.148 


> 0.7 


G3 


+2.052 ± 0.200 


2.397 ±0.152 


blue 


SDSS J1335+0118 




x[arcsec] 


y[arcsec] 


z 


A 


0.000 ± 0.001 


0.000 ± 0.001 


1.57 


B 


-1.038 ±0.002 


-1.165 ±0.002 


1.57 


G 


-0.769 ±0.011 


-0.757 ±0.011 


0.44 


SDSS J1353+1138 




R.A. (J2000.0) 


Decl. (J2000.0) 




A 


13 53 06.35 


+ 11 38 04.81 


1.63 


B 


13 53 06.08 


+ 11 38 01.43 


1.63 


G 


13 35 06.10 


+ 11 38 00.39 


~ 0.3 



Table 1. Astrometric properties and redshifts of the five lensed 
SDSS quasars: SDSS J0903+5028 from Johnston et al. (2003), SDSS 
J1001+5027 and SDSS J1206+4332 from Oguri et al. (2005), SDSS 
J1335+0118 from Oguri et al. (2004) and SDSS J1353+1138 from 
Inada et al. (2006). Units of R.A. are hours, minutes and seconds, and 
units of Dec. are degrees, arcminutes, and arcseconds. All x,y positions 
are relative to the A component and their errors do not include the error 
of the image scale. The positive direction of the those coordinates are 
to the west and north, respectively. 



3.4. SDSSJ1353+1138 

Inada et al. (2006) discovered the SDSS J1353+1 138 lensed 
quasar system from SDSS lensed quasar survey. For imaging and 
spectroscopy Inada et al. (2006) used the University of Hawaii 
2.2-m telescope, the Keck I and II telescopes and the Magellan 
Consortium's Landon Clay 6.5-m telescope. The observations 
showed that the two quasar components separated by l'/41 have 
redshift z = 1 .629 and the lensing galaxy in between the images 
is at z ~ 0.3. 



4. Photomery and image deconvolution 

4. 1 . Image deconvolution 

Time delay estimation requires high precision photometry sepa- 
rately for each lensing image. Since our targets are small angular 
separation systems we need to use a mathematical method, de- 
convolution, to separate the images. 

An observed image is a convolution of a real light distribu- 
tion with the so-called total blurring function or point spread 
function (PSF). The goal of the deconvolution is to extract 
the brightness of a source using knowledge about the PSF. 
Deconvolution is an inverse problem without a unique solution. 
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Fig. 2. Finding chart of the two systems: SDSS J1335+0118 and SDSS 
J1353+ 1138. The reference star and the star used to model the PSF are 
indicated 




SDSSJ1001+5Q27 




SDSSJ1206+433 




SDSS 11335 + 0118 



The method used for choosing the best solution is minimization 
of the difference between a model and a real image. 

All data used in this paper were deconvolved using the 
PSF Controled Deconvolution software created by Magain et al. 
(1998). In this software the input are the position and intensities 
of the images, which need to be well defined to succeed. We 
take the image positions from previous studies because these 
were better resolved and the quasar postions well determined 
(see Table 3.1). 

In Figure 3 we present the deconvolution results for the ob- 
served targets. On the left are shown the deconvolved images 
(mathematical models of lensing systems) and on the right the 
data images. We can see that the deconvolution reveals not only 
the quasar images but also the lensing galaxies which are not 
visible in the raw data. 



5DSS 11353+1138 



4.2. Photometry 

Photometry of all objects presented in this paper was made us- 
ing the PSF Controled Deconvolution software (Magain et al. 
1998). The software deconvolves all frames of a given object 
simultaneously, it constrains the position of the images and lens- 
ing galaxies and allows the magnitude of the quasar images 
to vary freely. The algorithm has been implemented in various 
analysis of quasars (Burud et al. 2000, 2002a,b; Hjorth et al. 



Fig. 3. Results of the deconvolution of the 5 systems. On the left are 
deconvolved images and on the right original images. The deconvolu- 
tion reveals not only quasar images but also the lensing galaxy which is 
otherwise not visible. 
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2002; Jakobsson et al. 2005; Vuissoz et al 2007; Eigenbrod et al. 
2006a). 

The light curves of SDSS J0903+5028, SDSS J1335+01 18, 
SDSS J1353+1138, SDSS J1001+5027, SDSS J1206+4332 
consist of 34, 31, 31, 24, 34 R-band data points, respectively, 
as presented in Figures 4-8. Data were taken from two optical 
cameras, ALFOSC (marked in black) and StanCam (marked in 
blue). The magnitudes are calculated relative to the PSF mag- 
nitude marked on the finding charts (Figures 1 and 2). Due to 
the small field of view and comparatively small sensitivity of 
StanCam and lack of bright stars next to our targets the refer- 
ence stars presented in the plots are only from ALFOSC. The 
reference star images were also deconvolved to look for system- 
atic errors in the deconvolution and to estimate the photometric 
uncertainties. Assuming that reference stars are not intrinsically 
variable any variability seen in their light curves must be due to 
photometric and deconvolution uncertainties. The error bars of 
the quasar images and the reference stars coming from ALFOSC 
are the averaged photon noise and uncertainties measured from 
the variations of the reference stars combined in quadrature. For 
data points from StanCam we have assumed one standard error 
bar of 0.02 mag which should approximately include all uncer- 
tainties. In the case of ALOSC data of SDSS J0903+5028 and 
StanCam data of SDSS J1335+01 18 there were no non-variable 
star in the field except PSF star, thus the error bars of the quasar 
magnitude are set to 0.02 mag., which is the approximate photo- 
metric uncertainty of the B image in average atmospheric condi- 
tions. 




100 150 
JD- 2454100 



200 250 



Fig. 5. R-band light curves of SDSS J1335+01 18. A (diamonds) and B 
(triangles) images of the quasar are plotted with black and blue colors 
indicating weather data coming from ALFOSC or StanCam. The ref- 
erence star is plotted with red squares. A polynomial is fitted to the A 
counterpart of the quasars light curve. This polynomial is also fitted to 
the B image with magnitude and time shift. The time shift is taken from 
the theoretical prediction of the time delay of this system (see Table 5). 



5. Light curves 

5.1. QSOs with little variability 

In 3 of the 5 systems (SDSS J0903+5028, SDSS J1335+0118, 
SDSS J1353+1138) no significant variability was observed. In 
Figures 4, 5 and 6 we plot the light curves of the three systems. 
By visual examination we see no variability in the quasar im- 
ages. 



B light curves with magnitude and time shift. Since there are no 
visible peaks in either light curves, the magnitude shift was cho- 
sen simply by taking a difference between average magnitude of 
image A and image B. The time shift is taken from the theoret- 
ical prediction of the time delay of those systems (see Section 
8). The fitted polynomials do not show any preferred time shift. 
This means that due to the lack of quasar variability during the 
time span of the monitoring we did not manage to measure time 
delays in these systems. 



20.2 
20.3 

I 20 4 

DC 

20.5 

20.6 
20.7 



A Image A 
c= Image B 



100 200 300 400 
JD - 2453600 



500 600 



Fig. 4. R-band light curves of SDSS J0903+5028. A polynomial is fitted 
to the light curve of A counterpart. This polynomial was also fitted to 
the B image with magnitude and time shift (see Table 2). 



In order to see whether small fluctuations in the light curves 
say something about the time delay we fitted polynomials to the 
light curves of the A counterpart. These were then fitted to the 



17.0 



17.5 



cn 
E, 

DC 



18.0 



18.5 




100 150 
JD - 2454100 



250 



Fig. 6. R-band light curves of SDSS J1353+1 138. Shapes and colors of 
the data point and fitting procedure identical with previous light curve. 
The time shift is taken from the theoretical prediction of the time delay 
of those systems (see Table 6). 
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5.2. Quasars with variability 

In 2 out of the 5 monitored gravitationally lensed quasars vari- 
ability was detected. In Figures 7 and 8 we show the light curves 
of SDSS J1001+5027 and SDSS J1206+4332, respectively. 



5.2.1. SDSS J1 001 +5027 

Figure 7 shows the light curves of SDSS J1001+5027 quasar 
images. We see that both images have some small variabilities 
in the light curve and also that both the A and B images have 
steadily decreased their brightness. 



17.5 



tS 
E, 
DC 



18.0 



18.5 



19.0 



A Image A 
c= Image B 

□ reference star - 1 ma 



ALFOSC 
STANCAM 




II 



50 



100 



150 200 
JD - 2454000 



250 



300 



A Image A 
o Image B 

□ Reference star +■ 0.5 ma 



19.0 



cn 
ts 
E, 

DC 



19.5 



20.0 - 



50 100 150 200 
JD - 2454050 



250 300 



Fig. 8. R-band light curves of SDSS J 1206+4332. Magnitude variabil- 
ities are seen in both light curves A (diamonds) and B (triangles) as 
bumps. The reference star is plotted with red squares. 



6. Time delay of SDSS J1 206+4332 

From the geometry of the system SDSS J 1206+4332 (see Figure 
3) we anticipate that the A image is leading since it is farther 
from the centre of the main lensing galaxy. Thus, we predict 
that any intrinsic quasar variabilities should appear first in the 
A image. From the mass modeling we know also that the time 
delay for this system is >50 days (see Section 8). We apply these 
constraints when measuring the time delay from the light curves. 



Fig. 7. R-band light curves of SDSS J1001+5027. Polynomials (solid, 
black lines) and linear regressions (dashed, red lines) are fitted to the A 
and B counterpart of the quasars separately. Shapes and colors of data 
point identical as described on previous plots. 



We have fitted a 7th order polynomial (see Figure 7 - black 
solid line) to the A data set. This was then fitted to the B light 
curve with magnitude shift which is the average magnitude dif- 
ference between light curves (0.4 mag) and time shift which is 
the predicted time delay (52 days). 

We also have fitted a straight line to light curves (see Figure 
7 - red dashed lines). The linear regression shows that both light 
curves have very similar evolution. Both images decreased their 
brightness by 0.2 magnitude during the first 200 days. This indi- 
cates that the brightness decrease is due to long intrinsic quasar 
variability (de Vries et al. 2003). 

Although variability is clearly visible, it is impossible to 
measure the time delay for this system just from the slope. The 
small fluctuations on the slope also do not give any conclusive 
results, neither visual shifting nor polynomial fitting help to find 
the time delay. 

5.2.2. SDSS J1 206+4332 

In Figure 8 we show the light curve of SDSS J1206+4332. We 
see clear long variabilities for both the A and B light curves. 
The variabilities are 100-150 days long, so the gaps in the sam- 
pling do not strongly influence the precision of the time delay 
estimation. The observed variabilities, consisting of bumps, al- 
low for measurement of the time delay for SDSS J1206+4332, 
as detailed in section 6. 



18.8 - 



€> lmageA + 116 shitted 
a Image B 



19.0 



CT) 

E, 

DC 



19.2 



19.4 



19.6 - 




100 200 300 

JD - 2454050 



400 



Fig. 9. R-band light curves of SDSS J 1206+4332. Magnitude variabil- 
ities are seen in both light curves A (diamonds) and B (triangles) as 
bumps. The A light curve is shifted in time by 115 days. The red line 
indicates the intervals used for cross correlation. 



The simplest form of such a measurement is visual shifting. 
Knowing that the variability has to appear first in the image A 
and then at least 50 days later in the image B we can associate 
the bump at the beginning of the A light curve centreed at ~50 
days with the bump in the B light curve centreed at ~150 days. 

A more quantitative measurement of a time delay is obtained 
by polynomial fitting. We fitted 7th order polynomials to both 
light curves using minimum^ 2 . In Figure 9 we show the two 
polynomials (solid lines). We have marked in red the parts of the 
polynomials which were used for calculating the time delay. In 
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CD 
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100 



110 120 
Time delay 
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Fig. 10. Results of cross correlation between data on one image and 
fitted polynomial fitted to light curve of second image. Solid lines rep- 
resents cross correlation between the A data and the polynomial fitted 
to B data; dashed lines are representing cross correlation between the B 
data and a polynomial fitted to the A data. 



order to calculate the time delay we have shifted the A polyno- 
mial in time and magnitude to fit it to the B data points and the B 
polynomial was shifted in time and magnitude to fit it to A data 
points. For each shift of the polynomials the goodness of fit was 
calculated. The results of these fits are presented in Figure 10. 
For magnitude shift, B - A — 0.3 1 mag the average of the two 
time delay estimations is 1 16^ days, where errors are uncertain- 
ties of the two values of the time delay combined in quadrature. 
Based on our modeling (Gl - SIE and G3 - SIS) (see Section 8) 
and using the measured time delay (At = 1 16^ days) we find a 

1 Mpc- 1 . 



Hubble constant of 73^j km s 



7. Microlensing 

Chang & Refsdal (1979) predicted that in lensed quasar systems 
the light path should be affected by stars in the lensing galaxy. 
Moving compact objects in the lensing galaxy can cause spec- 
tral changes, brightness variability and, in the case of multiple 
images, flux-ratio anomalies in the lensed quasar. 

The images of a lensed quasar may vary due to intrinsic 
quasar brightness changes and/or microlensing. Microlensing af- 
fects the light paths of each image differently (in the simplest 
case only one path is affected) whereas the intrinsic variations 
show up in all the images but at different times due to the time 
delay. Therefore, one can isolate the microlensing signal by sim- 
ply calculating the difference between two light curves. 

If we shift one of the light curves in magnitude (by the mag- 
nitude difference between light curves) and in time (by the time 
delay) and subtract it from the other quasar counterpart light 
curve the remaining variation in the light curve difference should 
not belong intrinsically to the quasar but rather to microlensing. 

Our previous study (Paraficz et al. 2006) showed that mi- 
crolensing variabiliy can be detected in this way. The monitoring 
program presented in this paper was designed to be sensitive to 
such microlensing signals. The only limitation are S/N and gaps 
in the light curve due to bad weather conditions. In Figures 1 1- 
15 we present the result of a microlensing variability analysis for 
our observed targets. 

The top plots present the quasar image lightcurves, with one 
of the images being shifted in time (by the theoretically pre- 
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Fig. 11. Microlensing of SDSS J1206+4332. Top: Time-delay shifted 
light curves, with the A image offset by 0.32 mag and 115 days. 
Bottom: Difference between linearly interpolated A image and B im- 
age. 
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Fig. 12. Microlensing of SDSS J1001+5027. Top: Time-delay shifted 
light curves, with the A image offset by 0.4 mag and 52 days. Bottom: 
Difference between linearly interpolated A image and B image. 



dieted or, if available, the measured time delay) and in magni- 
tude (by the average magnitude difference between the images). 
The lower panels show the light curve differences. The differ- 
ence is calculated between a linearly interpolated B light curve 
and the A data points. In Figures 11-15 we see that there is no 
signal stronger than 0. 1 mag which is ~ 2 ex. We therefore con- 
clude that we did not detect any significant microlensing event 
in any of the systems during our monitoring program. 



8. Modeling 

8. 1 . Analytical modeling 

In principle, one can calculate the Hubble constant from the time 
delay of a gravitationally lensed system, but there is one major 
obstacle, namely the mass-sheet degeneracy (Falco et al. 1985). 
This degeneracy is between the Hubble constant and the radial 
density profiles of the lens (Wucknitz & Refsdal 2001) and its 
angular structure (Freedman et al. 2001). Summarizing, differ- 
ent mass profiles give different results on the Hubble constant 
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Fig. 13. Microlensing of SDSS J1335+01 18. Top: Time-delay shifted 
light curves, with the A image offset by 1.35 mag and 26 days. Bottom: 
Difference between linearly interpolated A image and B image. 
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Fig. 14. Microlensing of SDSS J1353+1138. Top: Time-delay shifted 
light curves, with the A image offset by 1.1 mag and 17 days. Bottom: 
Difference between linearly interpolated A image and B image. 



without changing the lensing configuration. To break this de- 
generacy we need to have information about the mass profile, 
e.g., through the velocity dispersion (Romanowsky & Kochanek 
1999). For our systems we do not know which lens model is 
the correct one, but fortunately in many cases a lens can be 
quite well modeled using a fairly simple mass model like SIE 
(Singular Isothermal Ellipsoid) (Kormann et al. 1994) or NFW 
(Navarro-Frank- White) (Golse & Kneib 2002). 

In this paper, for the analytical modeling, we have modeled 
all our systems using singular isothermal potentials. We used 
those models because of their simplicity and because they agree 
with the physical properties of many observed lenses (Kochanek 
1993; Rusin et al. 2003; Koopmans et al. 2006). 

In the analytical modeling we have the positions of the two 
images as constraints and as free parameters we have the lens 
velocity dispersion, <t and the ellipticity, e. We also have two 
fixed parameters, the position angle, 6 e and the position of the 
lens galaxy. The number of degrees of freedom is 0, thus we can 
look for the one model which perfectly fits the data. The posi- 
tions of the images and lens are visible on the decovolved images 
(see Figure 3). The position angle of the lenses can be roughly 
estimated from analysis of the deconvoled images, as the mass 
profile of a lens tends to align to its visible component (Keeton 
et al. 1998). We allow this constraint to vary within 45°. We have 
not used flux ratios as constraints because of the possible influ- 
ence of reddening by dust (EKasdottir et al. 2006), microlensing 
(Paraficz et al. 2006) or small-scale structure in the lens potential 
(Dalai & Kochanek 2002). 

The ellipticity used here is defined as e = {a 2 - b 2 ) I {a 2 + b 2 ), 
where a and b are the major and minor axis. The position 
angle corresponds to the direction of the semi-major axis of 
the isopotential counted east of north. The modeling has been 
performed using the LENSTOOL software package available 
at http://www.oamp.fr/cosmology/lenstool/ (Julio et al. 2007). 
LENSTOOL is a software created for modelling strong lensing 
systems with parametric methods which chooses models using 
the Bayesian evidence. In all models we set the Hubble constant 
to be H = 73 km s" 1 Mpc" 1 . 

8.1.1. SDSS J0903+5028 
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Fig. 15. Microlensing of SDSS J0903+5028. Top: Time-delay shifted 
light curves, with the A image offset by 0.32 mag and 79 days. Bottom: 
Difference between linearly interpolated A image and B image. 



We have fitted the two simple models, SIS (Singular Isothermal 
Sphere) and SIE (Singular Isothermal Ellipsoid) to the lens- 
ing galaxy of SDSS J0903+5028. We have assumed that only 
one galaxy, the central one, influences the quasar light, even 
though other galaxies are visible to the north and south of the 
images (see Figure 16). The lensing images of the quasar are not 
collinear with the central galaxy which indicates that there is a 
quadruple moment in the potential. This moment can come from 
the tidal effect of the nearby galaxies or may be due to elon- 
gation of the central galaxy. Because of the degeneracy we are 
not able to distinguish which one of the possibilities is true, thus 
we use one lens, an SIE model which accurately predicts im- 
age positions. In Table 2 we summarize the results of the SDSS 
J0903+5028 modeling where the central single lens is an SIE. 
Based on this model, the predicted time delay for the system is 
79.4 days. The uncertainty in this value is entirely dominated by 
systematic modeling uncertainties. 



8.1.2. SDSS J1 001 +5027 

The lens of SDSS J1001+5027 consists of two galaxies (see 
Figure 3) which according Oguri et al. (2005) are at similar red- 
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Fig. 16. Lens model of SDSS J0903+5028. The figure shows an image 
of the system with over-plotted results of lens modeling. Blue and red 
lines represent critical and caustic lines, respectively. 
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Table 2. Lens model of SDSS J0903+5028. Modeled parameters are 
galaxy position, x,y, lens velocity dispersion, o"q, ellipticity, e, position 
angle, Q € and time delay, At. 



shifts 0.2 < z < 0.5. For our modeling we fix both lensing galax- 
ies at redshift 0.3. 

We have created two possible models of the lens environ- 
ment of SDSS J1001+5027. In the first one the system has only 
one galaxy with elongated SIE mass profile. This model is able 
to ray-trace the position of the quasar images. From Figure 17 
we also see that the position angle, P.A. of the model is consis- 
tent with the light distribution of the main galaxy but it is also 
in the direction of the other galaxy. This leads to the conclusion 
that the quadruple moment might come from the second galaxy 
rather then from the elongation of the main one. In the second 
version of the lens model we have set both gravitational lenses 
as Singular Isothermal Spheres, SIS. This model also reproduces 
the position of the quasar images and gives a sensible result of 
the galaxy velocity dispersions being 223 km s _I for main lens 
(Gl) and 159 km s _1 for G2. The results from both models are 
presented in Figure 17 and summarized in Table 3. The predicted 
time delay of the system is 34 days based on the SIE model and 
52 days based on the double SIS model. 
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Table 3. Lens model of SDSS J1001+5027. Parameters are as in 
Table 2. 
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Fig. 17. Lens model of SDSS J1001+5027. The figure shows an image 
of the system with over-plotted results of lens modeling. Blue and red 
lines represent critical and caustic lines, respectively. Top plot present 
the single SIE lens model and the bottom plot presents the two SIS lens 
models. 



8.1.3. SDSS J1 206+4332 

We have fitted SIS and SIE to the main lensing galaxy of 
the SDSS J1206+4332. We have assumed that the only galaxy 
which influences the quasar light is the central galaxy Gl (see 
Figure 3). Unfortunately we come to the same conclusions as 
(Oguri et al. 2005) that neither of these mass models gives physi- 
cally realistic results. The system seems to be more complicated, 
so presumably the other galaxies seen in vicinity of the system 
may play a role. A model of the main galaxy Gl (SIE) and G2 
(SIS) does not give any physical solutions either, which is ex- 
pected since the redshift of G2 is predicted to be rather small. 

It appears that the third galaxy lying ~ 3" north of the sys- 
tem (see Figure 3) has a strong influence on the system geom- 
etry. If we assume our model to be two lensing galaxies with 
one of them being Gl (SIE) and the other G3 (SIS) it is possi- 
ble to reproduce the positions of the images. The results of that 
final fitting are summarized in Table 4. Based on this model, 
the predicted time delay of the system is 1 12.5 days which is in 
agreement with our measurement of 1 16^ days (see section 6). 
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Fig. 18. Lens model of SDSS J 1206+4332. The figure shows an image 
of the system with over-plotted results of lens modeling. Blue and red 
lines represent critical and caustic lines respectively. 
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Table 4. Lens model of SDSS J1206+4332. Parameters are as in Table 2 
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Fig. 20. Mass model of SDSS J 1 353+ 1138. Figure shows an image of 
the system with over-plotted results of lens modeling. Blue and red lines 
represent critical and caustic lines, respectively. 
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Table 5. Lens model of SDSS J1335+01 18. Parameters are as in Table 2 



8.1.5. SDSS J1353+1 138 

SDSS J1353+1 138 has a centrally aligned lensing galaxy and 
the B image very close to the galaxy centre. There are no other 
close-by galaxies; thus the lens can be modeled with a simple 
SIE. The predicted time delay for that system is 17 days (see 
Table 6) which is in agreement with previous studies. 



# 


Type 


X 


y 


c 


<? e e 


At 






arcsec 


arcsec 


km s" 1 




days 


G 


SIE 


-0.25 


-1.04 


177 


0.46 35° 


17.2 



Table 6. Lens model of SDSS J 1353+1 138. Parameters are as in Table 2 



Fig. 19. Mass model of SDSS J1335+0118. Figure shows an image of 
the system with over-plotted results of lens modeling. Blue and red lines 
represent critical and caustic lines, respectively. 



8.1.4. SDSS J1335+0118 

The system has one lensing galaxy, at redshift z = 0.5, which 
we have modeled as an SIE. The lensing galaxy is barely vis- 
ible in our images, hence the position of the galaxy was taken 
from Oguri et al. (2004). According to observations made by 
Oguri et al. (2004) with Keck, the lensing galaxy is misaligned 
with the QSO images which indicates the existence of external 
shear or elongation of the lensing galaxy. A model with an SIE 
lens accurately predicts the image position and gives plausible 
galaxy parameters. The model details are summarized in Table 
5. The estimated time delay of the system is 26.2 days which is 
in agreement with previous studies. 



8.2. Non-parametric modeling 

We have also modeled our systems based on a non-parametric 
method proposed by Saha & Williams (1997) using the publicly 
available code PixeLens (Saha & Williams 2004). 

PixeLens reconstructs a pixelated mass map of the lens by 
implementing Bayesian statistics. It generates numerous lens 
models which fit the lens system geometry. There are two ma- 
jor advantages of using this method. First, it allows to ex- 
plore a wide range of models not restricted with parameters, 
which circumvents the non-uniqueness problem of lens mod- 
eling. Second, it provides an estimate of the systematic errors 
on the modeled time delay which is not possible to get when 
using analytical methods. Our purpose of using it is to check 
the robustness of our time delay estimations, i.e. we want to 
know whether the two approaches, analytical and pixeleted, give 
consistent results for the time delays. PixeLens generates mod- 
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Fig. 21. The modeling results of the five lensed quasars. On the left side 
we see potential contours modeled by PixeLens which are over-plotted 
on the mass map created in analytical modeling. On the right the are 
time delay histograms generated for each system by the 300 suitable 
lens models, arrows indicate the time delays estimated from analytical 
modeling. We see that the potentials of the PixeLens correspond quite 
well to the analytical models 



els which reproduce the exact position of lensing images and 
aligned within 45° of the visible lens P.A. All our systems were 
modeled assuming existence of external shear and all of them, 
except SDSS J1206+4332, were modeled as symmetric. For 
each system we generated 300 lens models which were used to 
estimate the average and median time delay and their model un- 
certainties. The modeling results of the 5 gravitationally lensed 
QSOs with PixeLens are summarized in Table 7. Figure 21 
shows the potentials and time delays of the 5 systems modeled 
with pixelated method. On the left side we see potential con- 
tours modeled by PixeLens which are over-plotted on the mass 
map created in analytical modeling. On the right are histograms 
of the time delays generated for each system by the 300 lens 
models. 

We see that the lens potentials of the PixeLens correspond 
to the analytical models and from the Table 7 we see that the 
estimated time delay values from analytical modeling lie within 
error bars of the PixeLens time delays. This is reassuring and 
gives confidence in our approach. We note also that the analyti- 
cally obtained time delays are systematically lower than the av- 
erage values obtained using PixeLens by 20-50%. 
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Table 7. Summary of modeled time delays of the five lensed quasars. 
Two different modeling methods are compared. For the pixelated 
method we generate 300 lens models for each system which are used 
to estimate the average and median time delay and their error bars. We 
see that the estimated time delays from analytical modeling are consis- 
tent with he PixeLens time delays. 



9. Simultanous modeling 

The Hubble constant estimated from time delay measurements 
is model dependent, as illustrated above. Different models give 
a different Hubble constant and measuring the time delay with 
high precision does not resolve that issue. One way of deal- 
ing with this problem is to generate a large number of mod- 
els of many lenses and create a distribution of possible val- 
ues for the Hubble constant based on these models. Thus, we 
use the pixelated method for several gravitationally lensed sys- 
tems constraining them to have the same shared value of Hq. 
Simultaneous modeling has been done in the past by e.g. Saha et 
al. (2006) - with 10 lenses and Coles (2008) - with 1 1 lenses. 

Figure 22 shows the result of the simultaneous modeling 
of 5 systems with time delays obtained at the NOT. These are 
RX J091 1+055 (Hjorth et al. 2002), SBS B1520+530, (Burud 
et al. 2002b), B 1600+434 (Burud et al. 2000) and J095 1+263 
(Jakobsson et al. 2005) and the time delay of SDSS J1206+4332 
reported in this paper. On the histogram the results of our simul- 
taneous modeling of the 5 QSOs are presented. For each object 
100 lens models were created. The Hubble constant estimated 
using this method is Hq = b\.5^ A km s Mpc -1 . 
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Fig. 22. The result of the simultaneous modeling, it is a histogram cre- 
ated from the modeling of the 5 QSOs monitored at the NOT for each 
of them 100 lens models were created, the Hubble constant estimated 
using this method is H = 61.5*4 km s Mpc -1 . 



10. Discussion 

We have estimated the time delay of the SDSS J1206+4332 to 
be At = 116^ based on a 280-day long monitoring campaign 
at the Nordic Optical Telescope. This shows the feasibility of 
measuring time delays from campaigns not much longer than 
the time delay itself. To model the system, we first assumed 
a single lensing galaxy modeled with the Singular Isothermal 
Ellipsoid (SIE), but were unable to reproduce the image posi- 
tions. By adding to the lens model the G2 and G3 galaxies we 
found that the influence of the G2 galaxy is negligible due to its 
small redshift but the G3 galaxy has a major impact on the sys- 
tem's geometry. A lens model with two galaxies, where the main 
one is a SIE and the second is a SIS, predicts the image positions 
very well. Thus this model resulted in a Hubble constant of 73^ 
km s Mpc~' assuming Q m = 0.3 and Q.^ = 0.7. 

Lack of short-term variability in the other monitored systems 
(SDSSJ 0903, SDSS J1001+5027, SDSS J1353+1138, SDSS 
J1335+01 18) meant that we did not succeed in measuring their 
time delays. From the quasar studies of Fohlmeister et al. (2008), 
Vuissoz et al (2007), Koopmans et al. (2000) and many others we 
see that during the 200-300 days of a quasar monitoring a quasar 
might be at the quiet stage or at the slow increase or decrease. 
Thus, despite our success in measuring the time delay of SDSS 
J 1206+4332, it appears that monitoring substantially longer than 
the predicted time delay of a given system will be unavoidable 
in many cases. 

We have modeled all five systems using both analytical and 
pixelated methods. Using the predicted time delays we have 
shown that none of the systems exhibited significant microlens- 
ing variability during the observing campaign. 

We have also performed simultaneous pixelated modeling 
with a common Hubble constant of five lenses for which time 
delays were measured at the NOT. The estimated Hubble con- 
stant from this analysis is Hq = 61.5*? km s _I Mpc~'. 
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Table.3. Photometry of two images of SDSS J0903+5028 quasar. 
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Q A/1 HI Q 

304 Hi 


20 


010052 


1 A "7 
10/ 


/2o29 


20 


306558 


20 


619706 


108 


77146 


20 


304743 


20 


618817 


109 


75729 


20 


303631 


20 


615568 


111 


"7 AAAA 

70990 


20 


302910 


20 


616263 


134 


49363 


20 


307745 


20 


615292 


135 


60198 


20 


305574 


20 


617903 


137 


56633 


20 


309453 


20 


618127 


138 


65267 


20 


305457 


20 


614339 


139 


65707 


20 


302961 


20 


618014 


157 


56273 


20 
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20 
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1/1 


c Anno 
54/ /O 


20 


305926 


20 
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172 
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20 
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OA 

20 
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20 


309605 


20 


6269/2 


191 


46353 


20 


308627 


20 


625028 


199 


52655 


20 


298319 


20 


617214 


428 


69561 


20 


1 A 1 O AO 

301808 


20 


627140 


430 


74676 


20 




20 


612383 


431 


69420 


20 


298253 


20 


6120/2 


432 


73686 


20 


299821 


20 


614314 


435 


74226 


20 


jUU /(J / 


20 


614775 


449 


59337 








(SI 


4S^ 


71216 


20 


296104 


20 


608069 


465 


70068 


20 


300696 


20 


615813 


467 


67382 


20 


298169 


20 


609918 


469 


74747 


20 


301014 


20 


613560 


483 


76553 


20 


300752 


20 


617619 


495 


56600 


20 


301731 


20 


615859 


501 


74698 


20 


303775 


20 


616146 


508 


74511 



Table.3. Photometry of two images of SDSS J1001+5027 quasar and a reference star. 



Rmag A 


Rmagg 


Rmag ref 


Err Rmag 


JD 


18.308300 


17.916424 


19.959659 


0.124381 


050 


69800 


18.323450 


17.937721 


19.929005 


0.034650 


062 


73200 


18.307941 


17.943148 


19.933155 


0.049381 


068 


69100 


18.277621 


17.921927 




0.022007 


074 


71047 


18.317622 


17.978474 


19.906741 


0.046542 


100 


66100 


18.340800 


17.996533 


19.931855 


0.028322 


103 


73100 


18.343249 


18.001368 


19.906263 


0.018813 


105 


76800 


18.342268 


18.016725 


20.028774 


0.028204 


109 


72400 


18.353127 


18.013164 


19.943019 


0.147739 


111 


65300 


18.346979 


17.990860 


19.911710 


0.010532 


123 


78600 


18.339169 


17.959495 




0.034687 


134 


75378 


18.387632 


18.042124 


19.966310 


0.011901 


138 


58600 


18.385108 


18.037049 


19.931430 


0.012798 


142 


76500 


18.415253 


18.066928 


19.936043 


0.018701 


155 


43300 


18.350806 


17.993212 




0.027625 


156 


58644 


18.363621 


18.000546 




0.020099 


158 


59506 


18.360164 


17.996483 




0.010037 


160 


60051 


18.409971 


18.072709 


19.948912 


0.013566 


176 


41200 


18.394950 


18.070155 


19.918187 


0.010659 


182 


41100 


18.435319 


18.079750 




0.037845 


215 


53896 


18.421549 


18.055671 


19.896153 


0.010842 


219 


38300 


18.417204 


18.075555 




0.037845 


221 


36005 


18.460061 


18.094696 


19.930324 


0.038369 


227 


41000 


18.459076 


18.113779 


19.938801 


0.093712 


228 


37800 


18.405639 


18.072108 




0.063487 


235 


38689 


18.409954 


18.081878 




0.018256 


239 


39029 


18.459080 


18.126627 


19.937951 


0.010008 


247 


42200 


18.440717 


18.107298 


19.904796 


0.010067 


254 


39000 


18.429906 


18.072435 




0.091177 


260 


41507 


18.487492 


18.093200 




0.066955 


276 


40500 



Danuta Paraficz: Results of optical monitoring of 5 SDSS double QSOs with the Nordic Optical Telescope 
Table.3. Photometry of two images of SDSS J1206+4332 quasar and a reference star. 



RmagA 


Rmag B 


Rmag re f 


Err Rm3.g 
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18 


967233 


1 Q A G7 1 C ^ 

iy.4y / loj 


n ni n^o i 
U.UlUJol 


ncn AA7HA 
UjU.OO /UO 


iy.4jzoyj 


18 


985873 
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ly 
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19 


152/32 


iy.jo /y4o 


n nQ7Q7^ 
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iy.444USz 


19 


135932 
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ly.jZz III 
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19 
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iy.4/o4y / 


19 
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19 
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m /i A7A7n 
iy.4o/u /u 
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iy.4cuuuy 


19 


133444 


iy.jzo4Uj 
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ZZ0.41 I/O 


1 Q A OQ 1 QO 

iy.4oy 1 jo 


19 


153896 


1 Q ^AA 
iy.j44?53o 


U.U4j j 1 1 


OQA /I 1 777 

ZjU.41 / / / 
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iy.4yyj4U 


1Q 
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151 204 


iy.3 / /4y2 


U.1Z3JJ / 


m /ii i/ic 
zjz. 41 145 


19.475140 


19 


120366 


19 341953 


159005 


234 46066 


19.490768 


19 


109479 


19.534108 


0.034926 


240.41984 


19.506084 


19 


096243 


19.498606 


0.010205 


244.40898 


19.500106 


19 


028428 


19.390230 


0.110865 


265.40616 


19.449340 


18 


985593 


19.512254 


0.015322 


269.39189 


19.502726 


19 


039390 


19.473984 


0.028473 


281.37866 


19.507411 


19 


054861 


19.307993 


0.192910 


283.37888 



Table.3. Photometry of two images of SDSS J1335+1 138 quasar and a reference star. 



RmagA 


Rmage 


Rmagref 


Err Rmag 


JD 


18 


103033 


19.366769 




0.020000 


30.724810 


18 


173849 


19.445618 




0.020000 


35.773720 


18 


168048 


19.507642 


19.811958 


0.017249 


39.739070 


18 


138333 


19.480296 


19.827976 


0.010191 


42.774090 


18 


158132 


19.434327 




0.020000 


43.732740 


18 


116399 


19.338482 




0.020000 


49.774310 


18 


128543 


19.508473 


19.835769 


0.013971 


55.622420 


18 


132991 


19.461828 




0.020000 


61.534550 


18 


094535 


19.430734 




0.020000 


67.730140 


18 


109781 


19.461155 


19.816665 


0.013688 


100.67228 


18 


105817 


19.439126 




0.020000 


106.51100 


18 


109180 


19.381934 




0.020000 


110.62587 


18 


098207 


19.397216 




0.020000 


114.58660 


18 


110809 


19.398308 




0.020000 


116.54903 


18 


137283 


19.435018 




0.020000 


120.53586 


18 


099405 


19.501067 


19.810479 


0.018473 


126.52229 


18 


103190 


19.397234 




0.020000 


134.55204 


18 


086870 


19.403923 




0.020000 


136.39267 


18 


101415 


19.440893 




0.020000 


138.40052 


18 


094541 


19.434533 




0.020000 


140.38453 


18 


095505 


19.473765 


19.828225 


0.010241 


155.41550 


18 


101533 


19.475559 


19.835056 


0.013482 


157.45223 


18 


076486 


19.352538 




0.020000 


164.41317 


18 


082410 


19.449831 


19.824963 


0.010035 


179.41976 


18 


086675 


19.448914 


19.826855 


0.017975 


181.40853 


18 


088464 


19.443357 


19.840949 


0.011272 


183.41002 


18 


063292 


19.420842 


19.820810 


0.015805 


199.42655 


18 


017446 


19.390862 


19.838252 


0.080799 


213.41659 


18 


052497 


19.417072 


19.745833 


0.012022 


214.39304 


18 


045362 


19.375354 


19.832686 


0.020918 


228.41725 


18 


069432 


19.327107 


19.844385 


0.016030 


232.38603 


18 


057617 


19.398826 


19.838541 


0.038127 


234.37637 



Danuta Paraficz: Results of optical monitoring 



of 5 SDSS double QSOs with the Nordic Optical Telescope 



Table.3. Photometry of two images of SDSS J1353+01 18 quasar and a reference star. 



RmagA 
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Rmag ref 
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